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Abstract 

The manipulation of trapped charged particles by electric fields is an accurate, robust and reliable 
technique for many applications or experiments in high-precision spectroscopy. The transfer of the 
ion sample between multiple traps allows the use of a tailored environment in quantum information, 
cold chemistry, or frequency metrology experiments. In this article, we experimentally study the 
transport of ion clouds of up to 50 000 ions. The design of the trap makes ions very sensitive to 
any mismatch between the assumed electric potential and the actual local one. Nevertheless, we 
show that being fast (100 ps to transfer over more than 20 mm) increases the transport efficiency 
to values higher than 90 %, even with a large number of ions. For clouds of less than 2000 ions, a 
100 % transfer efficiency is observed. 
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I. INTRODUCTION 


An ion cloud confined in a linear radiofrequency (RF) quadrupole trap is an example of 
a non-neutral plasma (NNP), a plasma consisting exclusively of particles with a single sign 
of charge. Its thermal equilibrium state has been studied in detail by Dubin, O’Neil and 
coworkers in the context of large ensembles in Penning traps (for a complete review, see 
[1]) and extrapolated to ions in RF-quadrupole traps [2, 3] and even to ions in multipole 
traps [4], Thanks to Doppler laser cooling, the transition of an ion cloud to a liquid and a 
Coulomb crystal phase was observed [5, 6] after it was theoretically studied [2, 7]. 

In this paper, we address an out-of-equilibrium issue with the transport of ion clouds 
by the translation and deformation of the trapping potential. Our aim is to shuttle ions 
between separate trapping zones, and the objective is to do this as fast as possible without 
loss. This problem is related to the transport of single ions in micro-traps, which has been 
realized mainly for quantum information processing (QIP) applications. In our experiment, 
a more general problem is studied as the transported ensemble is a many-body system with 
long-range interaction. However, we can compare the center-of-mass motion of the ion cloud 
to the trajectory of a single ion, and we therefore use the QIP transports as model systems. 

Different groups have addressed the question of transport of single ions [8-10] as it is 
a crucial issue for scalable architectures of QIP in ion traps. One of the main concerns in 
these experiments is to avoid heating issues during the transport, the implementation of 
robust gate operations also requires very high transfer efficiencies [11]. In micro-traps, the 
transport distances for a single ion are of the order of lOO/rm. Care is taken to translate 
the ion in a quasi-constant potential well, which requires a large number of electrodes to 
tailor the trapping potential along the transport path. Speed is an additional issue which 
has to be taken into account, as shuttling ions between different sites is only a preparatory 
or intermediate task in a more sophisticated protocol and should not last longer than the 
computational gate. 

Our experiment comprises two macroscopic linear quadrupole RF traps, storing ion clouds 
of a thousand up to a million ions. The trapping and shuttling along the common trap ; 2 - 
axis are controlled by three DC-electrodes. Such large traps are typically used in frequency 
metrology in the microwave domain [12], exotic ion studies [13] or experiments in physical 
chemistry [14]. Two different trapping zones are useful to keep one zone free from contact 
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potential induced by neutral atom deposit or to accumulate ions in one of the trapping zones. 
High efficiencies for transfer are mandatory, and rapid transport protocols can reduce dead 
times which are detrimental to frequency stability in the case of atomic clocks and to sample 
conservation in case of short-lived species. 

Many-body transport is also a concern for experiments with an ensemble of cold neutral 
atoms and Bose-Einstein condensates (BECs) which have been transported without heating, 
making use of shortcuts to adiabacity [15, 16]. In [15] a cloud of a few 10® cold atoms is 
shuttled back and forth with an optical tweezer over a distance of 22.5 mm, in times as short 
as four trap oscillation periods. The use of an optical tweezer is very advantageous as it can 
be moved without deformation and the faster than adiabatic transport scheme relies on this 
non-deformation. The scheme designed in [16] allowed the authors to translate a cold gas in 
the non-interacting limit as well as a BEC and the non-condensed fraction by more than half 
a millimeter. The transport and decompression of the atomic sample was engineered using 
dynamics invariants. Because of the Coulomb repulsion, the method used for cold atoms can 
not be extrapolated to ion clouds. For most experiments, adiabatic transfer is not a relevant 
solution. For the trapping potentials in our experiment, the adiabatic transfer time is of 
the order of several tens of seconds, incompatible with a majority of precision experimental 
protocols. 

This article is organised as follows : the experimental set-up and techniques are presented 
in section II. Section III is devoted to results and analysis of the transfer efficiency and the 
heating induced by the transport is analysed in section IV. Section V deals with the ion 
number effect for a specific transport duration with a focus on smaller clouds. Conclusion 
of this work can be found in section VI. 


II. DESCRIPTION OF THE EXPERIMENTAL SET-UP AND PROTOCOL 

A. Trapping and laser cooling 

Calcium ions are trapped in a two-part linear quadrupole trap of inner radius tq = 
3.93 mm which is designed to have reduced non-harmonic components in the trapping po¬ 
tential : the RF-electrodes are four cylindrical rods of total length 58 mm connected in 
a balanced way (±(VRir/2) cos(r2t)) to the RE supply (no grounded electrodes). The RF- 
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frequency, f2/27r, is 5.23 MHz and the potential difference between neighbouring electrodes 
oscillates with a peak-to-peak amplitude of 1045 V (if not mentioned otherwise) which gives 
a Mathieu parameter of = 0.15, well within the adiabatic approximation regime where 
the RF-trapping can be approximated by a static harmonic potential (the pseudo-potential) 
with frequency = gxf^/(2\/2) = 27r x 277 kHz. 

Trapping along the symmetry axis is reached by DC-voltages applied to electrodes per¬ 
pendicular to the rods. Three DC-electrodes are located at equal distances along the rods, 
creating two distinct trapping zones of 21 mm length, each. This double well conhguration 
is used for accumulation in one of the wells before further transport of the ions to other 
traps, in line with the quadrupole one. In [17] our protocol for genuine accumulation of 
ions is described in detail. In order to laser-cool the ions in both trapping regions using 
the same laser beam, the DC-electrodes must leave the trap x-axis free, which justihes the 
open shape of the three of them. Their design, detailed in [18], results from a compromise 
between reduction of the non-harmonic contributions in the potential and of the screening 
effect induced by the RF-rods. 

Calcium ions are produced by photo-ionization of neutral calcium atoms from an effusive 
beam crossing the trap axis perpendicularly in the horizontal plane. The photo-ionization 
process implies two photons and the hrst step (423 nm) is a resonant excitation tuned to 
select the most abundant isotope, ^°Ca [19, 20]. The second photon, at 375 nm, takes the 
atomic system above the ionization threshold. Both beams co-propagate along the trap axis. 
Ions are laser-cooled by two collimated 397 nm-beams on the [4S'i/2 — 4Pi/2]-transition, of 
equal power (2 mW on a 2 mm 1/e^ diameter), counter-propagating along the trap axis. 
Once excited from the ground state, calcium ions can relax to a long-lived metastable state 
[ 3 ZI 3 / 2 ] from which they have to be re-pumped to maintain efficient laser cooling. This re¬ 
pumping process is assured by a 866 nm laser beam [ 3 D 3/2 — 4 P 1 / 2 ] of approximately 2.5 mW 
and 4 mm 1/e^ diameter which co-propagates with one of the cooling lasers. Simultaneous 
ion creation and cooling allow to trap clouds as large a several hundreds of thousands of ions. 
For the work presented here, we tuned the ion creation parameters to reach a cloud size of 
the order of 20 000 ions, which takes typically an integration time of 15 seconds for laser 
powers of 80 /iW at 423 nm and 4 mW at 375 nm. Ions are detected by their laser-induced 
fluorescence at 397 nm which is collected through a dedicated objective with anti-reflection 
coating and a high numerical aperture (Sill Optics, / = 66.8 mm, N.A = 0.28). 
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FIG. 1. (a): Scheme of the experimental set-up and (b): calculated DC- potential along the trap 
axis if lOOOV are applied to each electrode. The distance between electrode centers is 23 mm and 
the thickness of each electrode is 2 mm. 


The distance L between the center of adjacent DC-electrodes is 23 mm and the trapping 
along the z-axis resnlts from the electric potential gradient of the snm of each electrode 
contribution. A solution of the Laplace equation by a finite difference method software 
(Simion [21]) gives the potential profile associated to each electrode of thickness 2 mm. 
The characteristic shape of the resulting potential wells along the trap 2 ;-axis are shown in 
Figure 1. The effective potential well, which can be estimated by the difference between 
the maximum and minimum total contribution, is lower than the DC voltage applied to the 
electrodes because of the screening by the RF-electrodes [18] and by the overlapping of the 
potential prohles at the center of the trapping zones which offsets the potential minimum. 
The large distance between electrodes results in a very small overlap of the profiles and as a 
consequence, in a low potential minimum. Nevertheless, the screening effect sets a limit to 
the axial trapping efficiency : for 1000 V applied to each electrode, the voltage calculated 
at the electrode center is 29 V for the central one and 28.8 V for the end electrodes whereas 
the minimal potential value is 1.16 V. The small overlap which leads to a deep potential 
well is a drawback for transport issues, as explained in the following. 
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B. Transport protocol 


As pointed out in the introduction, single ion transport without heating is a major issue 
for the scalability of trapped ion based quantum computer and is the subject of several 
experiments [8-10]. In these works, the authors characterize and compare different transport 
protocols with respect to the heating they induce on the ion motion. To guide us with 
our experiments, we have extrapolated these analysis for an ion cloud by using molecular 
dynamics (MD) simulations in [22]. We use the notations of previous work and we call 
(f)i{x,y,z) the electric potential created by the DC-electrode i when IV is applied to it. 
Then, the total DC potential inside the trap can be expressed as [23] 

N 

$(t, X, y,z) = Y^ Vi{t)(j)i{x, y, z) (1) 

i 

if V(t) is applied to electrode i. 

As the laser-cooled ions explore less than a tenth of the radial trap extension, we assume 
that the dependence of the DC-potential with the coordinates x and ?/, perpendicular to the 
trap axis, is not relevant and we call (pi^z) the on-axis evaluated function cpi^x = 0 , ?/ = 0 , 2 :). 
With our electrode geometry, the electric potential spatial distribution created on the axis 
is very well htted around its maximum by the equation : 

/.w = o. (1 + ‘ (2) 

with Zi the position of the center of electrode-i, Wi = 8.9 mm and 01,3 = 28.8 mV and 
02 = 29 mV when 1 V is applied on the i-electrode. The non-symmetric environment of 
the trap explains the small variation between the a* values. If \z — Zi\ = Wi, fi{z) = Oj/lb 
and we can consider that 2wi is a good enough approximation of the effective width of the 
potential prohle. When the three DC-electrodes are connected to the same potential, the 
potential well can be considered as harmonic around each potential minimum, behaving like 
mu‘^{z — ZcY/2. The value of Uz deduced from a £t of the potential around its minimum 
is 2ti X 107.5 kHz for Vj = 1000 V whereas a measurement by parametric excitation of 
the ion cloud [24] gives 27r x 124 kHz. The difference between these two values shows that 
the potential deduced from the calculation does not reproduce exactly the real potential 
experienced by the ions. 


6 



The transport protocol relies on the time variation of the Vi potential, designed to make 
the potential minimum obey a time prohle Zmin{t). This condition translates into 


(9$ 

dz 


.it) 


0 


(3) 


There are two local minima which meet each other when they reach the centre of the cen¬ 
tral electrode and the challenge of the experiment is to design a potential evolution which 
transfers the ions from the minimum in zone 1 to the minimum in zone 2 (see hgure 1). In 
the following, we call Zminif) the path we want the ions to follow and it can be written like 
[ 10 ] 

Zmin{t) = g{t){H{t) - H{t - tg)) + LH{t - tg) - Lj^ (4) 


with L the shuttling distance, g{t) the time prohle of the transport, tg its duration, and 
H{t) the Heaviside step function. Guided by numerical results detailed in [22], we used for 
the experiments presented here the time prohle described by 

L /tanh(4(2t/t3 - 1)) 

= 2 (-- 



These simulations showed that among four compared time prohles, the one following Eq. 5 
is the most robust against transfer duration variations. They also give evidence that the 
deformation of the trapping potential along the transport is responsible for the heating of 
the center of mass motion and of the ions’ motion in the center of mass frame [22] . In our 
experimental set-up, the distance between electrodes is larger than the ehective width of the 
potential they create, therefore, keeping the axial potential undeformed while translating 
its minimum requires huge voltages that we cannot provide. With the applied voltages Vi 
limited to 2000 V, we can only change the depth of the effective harmonic potential but 
cannot compensate for its deformation. 

In practice, the potential minimum is forced to obey the time prohle Zmin{t) if 


V2{t) 


Vi{t)4>'i{z) + Vz{t)(t)'^{z) 

4 ( 2 :) 


^min 


it) 


( 6 ) 


and the harmonic contribution of the resulting axial potential along the transport can be 
deduced by 


= 




m dz"^ 


.it) 


(7) 


where Q is the ion charge and m its mass. 
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Computing V 2 (t) through Eq. 6 requires to know the potential prohle created by each 
electrode to estimate the hrst order derivative Single-ion experiments have shown the 

great sensitivity of the transport induced heating on the precise knowledge of the potential 
geometry. Reference [25] proposes a characterisation method for micro-traps where a single 
ion explores the potential and its flaws. The size of our trap and of the ion sample are 
not suited for this method, and we have hence based our calculations on (pi^z) ~ fi{z) (see 
Eq 2). Eq ( 6 ) leads to a discontinuity of V 2 (t) for z^inif) = Z 2 i the center of the electrode 2 
where 02 (^ 2 ) = 0. To avoid this discontinuity, a constant relation between Viit) and V^{t) 
is imposed, given by 

U(i) = -U(i)0k4 <*> 

In our experiments, Vi and V 3 are kept constant. In a perfectly symmetric device Vslt) = 
Vi{t) would solve the problem but any asymmetry in the electrode environment breaks this 
equality. The reader is referred to [22] for details about how to avoid discontinuities in 
numerical simulations. 

III. TRANSFER EFFICIENCY FOR LARGE ION CLOUDS 

A. Estimation of the number of ions 

We are primarily interested in the transport efficiency as a matter of the relative number 
of ions passing from trapping zone 1 to zone 2. A precise quantitative study requires the 
measurement of the number of trapped ions, which is of the order of a few tens of thousands. 
The ion’s fluorescence signal is split between a photomultiplier and an intensified CCD 
camera. For hxed laser frequencies and trapping parameters, the fluorescence counting 
rate depends on the number of ions in the trap but also on their temperature. Because 
the transport can induce heating and ion loss, and because RF-heating depends on the 
ion number and their temperature [26], there is no simple relation between the recorded 
fluorescence signal and the number of trapped ions. Indeed, we very often observe a signal 
increase when the ion number has decreased, because of a smaller RF-heating. 

To develop a quantitative diagnostic independent of the signal counting rate, we use the 
density characteristics of the liquid phase of an ion cloud. One can show that in the cold fluid 
limit, a singly-charged sample in a harmonic potential has a uniform density [27], bound by 



an ellipsoid of revolution where the density falls to zero on the scale of the Debye length. 
This results from the Boltzmann-Poisson equation in the low temperature limit and for ions 
in a linear quadrupole trap the density depends only on the trapping pseudo-potential [4], 
This property, as well as the predicted aspect ratio of the ellipsoid [28], have been verihed 
quantitatively very accurately for ions in a linear quadrupole trap in [29]. 

Every ion ensemble is cooled to the liquid phase before and after transport in order to 
quantify the ion number. The difference between liquid and gas phase is easily detected 
by the variation in the fluorescence level [29]. The calculated density in the liquid phase is 
1.40 X 10^ mm“^ and for a typical temperature of 100 mK, the Debye length is 1.85 pm, 
which hts within two pixels on the camera with an optical magnihcation of the order of 
13. The typical size for the ion cloud in the liquid phase is of the order of 700 pm for the 
semi-major axis and 300 pm for the semi-minor axis of the ellipse. The Debye length is 
then negligible compared to the cloud size and we consider a uniform density all over the 
cloud in the liquid phase. By changing the laser cooling efficiency and thereby the sample 
temperature, we checked that, for ion clouds smaller than 50 000 ions, once the ellipsoid 
is formed by laser cooling, the measured dimensions for the ellipse are independent of the 
fluorescence signal level and the slight modihcation of the Debye length has no impact on 
the measured values. For larger clouds, transition to the liquid phase requires to reduce the 
RF-voltage amplitude and the border of the ellipses spread over a larger scale. We have 
developed a software which automatically hts and extracts the dimensions of an ellipse from 
the recorded picture of the collected fluorescence and the method used by this software is 
now developed. 

The first step is to dehne the contour of the fluorescence signal. This requires the dehni- 
tion of a threshold for the signal, independent from the number of photons scattered per ion. 
The threshold criteria is provided by the analysis of the section of the signal along one pixel 
line across the image. The derivative of the signal with respect to the pixel position shows 
two sharp extrema, Xi and X 2 at the ellipse border (see hgure 2). Their position falls in 
the same pixel as the one chosen by a £t ”by the eyes” and does not depend on the absolute 
level of the signal. As some of the ellipse pictures are longer than the camera detector, we 
base our protocol on a section located approximately on the small axis and define the level 
threshold as the mean signal (5'(Xi) -|- S{X 2 ))/ 2 . 

The second step is to fit the contour by an ellipse equation, including a possible angle 
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FIG. 2. (a) example of the picture of an ion cloud (dimensions in pixels) in the liquid phase with 
the fit of the contour (exposure time 0.5 s). The trap axis and laser propagation direction are 
vertical, (b) Smoothed section of the signal along the short axis of the ellipse, and (c) its derivative 


between the ellipse semi-axis and the detector rows. This is done through a 2D £t subroutine 
and produces a £t which falls in the same pixels as the original contour, which conhrms the 
ellipsoid shape of the cloud. To valid our £t procedure, for the same ion cloud, simply 
deformed by changing the trapping parameters, we use the two calculated semi-axis lengths 
Re and Le to estimate the ellipse aspect ratio pe = R^Le and volume 14 = AnRlLe/S 
extracted from experimental data. We can compare the aspect ratio with the expected one, 
deduced from the effective pseudo-potential p = fipjz/oJr) by an equation demonstrated in 
[28] and experimentally confirmed in [30]. More precisely, we measure the relative shift 
between the length deduced from the fit and the one deduced by the htted radius, assuming 
a known aspect ratio : 5l = {Le — Re/p)/Le. Furthermore, we can check that the volume 
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modifications obey what is expected from the density n(r) at low temperature : 

n(r) = eoAd'pp(r)/Q (9) 

where thpp(r) is the harmonic pseudo-potential, characterised by 

'J'pp(r) = + y"^) + (10) 

which leads to a uniform density n(r) = uq = eo2mul/Q‘^. As we want to reach a 1% 
level accuracy in relative volume estimation, we have to go beyond the hrst order adiabatic 
approximation. By expanding the calculation of the coefficient in the Mathieu solutions to 
the second order in (g^, a^), one can show that 

■ ( 11 ) 

In our case, Qx is induced by the 2 :-axis trapping-voltage Vdc, cix = and with 

our operating parameters, the correction is in the 1% range. Taking that into account, for 
the same RF amplitude but different DC-voltages, we observe relative fluctuations lower 
than ±1% for both the length and the volume of the same cloud in zone 2. In zone 1, 
we observe volume fluctuations that can reach 6 %, far larger than the length fluctuations 
6l which remain in the ±1% range. We attribute this difference to the contact potential, 
identified in zone 1, and induced by calcium deposition on the quadrupole rods, in front 
of the calcium oven. For different trapping parameters, the cloud is displaced in the trap, 
giving an optical image with a slightly different size. This is in particular true when the 
RF amplitude is changed and where apparent ion numbers can vary by 10%. For constant 
trapping parameters, like used for estimating the relative number of ions after a transport 
protocol, the uncertainty on the volume is ±1% for zone 2 and ±1.5% for zone 1 {6a 
conhdence). Precise investigations of transport efficiencies in terms of ratios of ion numbers, 
only requires measurement of the volume of the ellipses as we compare cold clouds with 
identical density. As for an estimation of the number of ions in a cloud, the fluctuations of 
the apparent particle number for the same cloud, when it is deformed and shifted, lead us 
to £x a 5% uncertainty on the absolute number (and the uncertainty induced by the optical 
magnification is negligible here). 
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B. From one local minimum to the other 


In a first step, we study our ability to transfer ions from one trapping zone to the other, 
depending on the duration of the transfer tg for given trapping parameters. Bandpass 
limitations of the DC-supplies prevent us from investigating transport durations shorter 
than 80 /xs. All the experiments were done with an ion cloud with a typical size ranging 
from 5000 to 20000 ions. For some well identihed transport protocols, we checked that the 
transport efficiency is independent of the ion number as long as this number is larger than 
2000. For clouds smaller than 2000, the efficiency is higher than for larger clouds and can 
reach 100 %. A focus on smaller cloud transfer is presented in section V. As our detection 
is based on the observation of induced fluorescence, the cooling laser beams remain applied 
during the transport. For some specihc transport protocols, we compared the transport 
efficiency with and without the cooling laser during transport. The observed differences 
were only of the order of a few % showing that the cooling effect does not play an important 
role. Indeed, the capture range of the Doppler laser cooling is 9.2 m/s, smaller or far smaller 
than the average shuttling velocity which ranges from 20 to 200 m/s. 

Figure 3 shows the fraction of ions leaving the trapping zone 1 for tg between 80/is and 
2.6 ms. The hrst major observation is that the number of leaving ions depends strongly 
on tg, alternating between nearly 0 and 100% several times before these oscillations are 
damped. Changing the axial trapping potential by changing the DC voltages, shifts these 
oscillations with tg and makes more oscillations observable for a larger value of the DC 
voltages. If the same protocol is applied to a smaller ion cloud (typically 1000 ions and 
smaller), oscillations are also observed with identical temporal imprint, excluding a number 
dependent effect. A possible explanation for this interchange between a high and low transfer 
probability is the oscillation of the ions from zone 1 to zone 2 and their return to trapping 
zone 1 before the transport function is completed. This assumption is tested by setting 1% to 
-2000 V during the transport protocol, which means opening the second trapping zone. For 
this conhguration, oscillations are still observed, shifted in tg with respect to the previous 
symmetric conhguration. These observations are consistent with the hypothesis that ions 
that were still in trapping zone 1 once the transport protocol was completed, entered zone 2 
far enough to see their trajectories modihed by the deformation of the potential, but not far 
enough to be attracted by the negatively polarised electrode. For transport durations longer 
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FIG. 3. Fraction of ions leaving the trapping zone 1 vs the transport duration tg for different 
on-axis DC-voltages : (a) Vi = 600 V, (b) Vi = 1200 V. The lines are a means to guide the eye. 


than 1 ms, the fraction of transferred ions reaches a stationary valne. In an ideal, symmetric 
system this valne is expected to be 50 %. In onr dnal trap, an asymmetry, very probably 
dne to the contact potential in trapping zone 1, can be responsible for this imbalance in ion 
repartition for long transport. 

To get more insight into this issne, we nse a MD simnlation to compnte the trajectory 
of a single ion in the translated and deformed potential applied in the experiment, as a 
good approximation of the center of mass motion [22]. Actnally, the resnlts of this simnla¬ 
tion depend on the method nsed to describe the potential. If the eqnations of motion are 
integrated in the calcnlated axial potential htted by Eq. 2, the probability for the ion to 
be transferred to zone 2 is nnity, whatever is the transfer dnration. To come closer to the 
experimental sitnation, we keep the waveforms Vj(f) as nsed in the experiments bnt inte¬ 
grate the eqnation of motion in the potential grid calcnlated by the hnite difference method 
software, Simion [21]. In this condition, oscillations of the probability to transfer the ion 
to trapping zone 2 are observed, depending on the transfer dnration. The conseqnence of 
the discrepancy between the two descriptions of the potential is visible on fignre 4 where 
several examples of ion trajectories are ploted. For tg = 189 ps, the ion trajectory follows 
the potential minimnm from trapping zone 1 to zone 2. For longer transfer dnrations, the 
ion is ahead of the potential minimnm and for tg = 884 /is, makes a U-tnrn in zone 2 before 
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FIG. 4. Single ion trajectories computed by MD, versus the relative time elapsed during the 
transfer. The trajectories are computed in the potential grid provided by the Simion software, 
based on our electrode geometry, when V 2 (t) obeys Eq. 6 computed with (l)i{z) = fi{z). Red 
dashed line : tg = 189 /rs, solid blue line : tg = 884 /rs, green dot-dashed line : tg = 1621 /rs. 
The black dotted lines are the position of the two potential minima along the transfer. Horizontal 
dash-dotted lines indicate the positions of the centres of the traps. 

ending in zone 1, like assumed previously. For even longer transfer times {tg = 1621 ps), the 
calculated trajectory shows a U-turn in zone 2, followed by a U-turn in zone 1 to hnally have 
an ion efficiently transferred to zone 2. Longer transfer durations give rise to an increasing 
number of U-turns which results in an oscillation between trap 1 and trap 2 for the ion hnal 
position. The experimental results exhibit a shorter timescale than the simulations. This 
can be explained by a larger difference between the potential assumed to compute the Vi{t) 
and the one experienced by the ions in the trap. The inconsistency identihed above cannot 
be avoided in a large-scale experiment where it seems unrealistic to generate a precise map 
of the complete electric held seen by the ion cloud along its transport. Nevertheless, the 
experimental results show that it is still possible to force the ion to transfer even if the cor¬ 
responding time scale can not be exactly foreseen. The next step in the transport efficiency 
analysis is to look at how many ions effectively settle in the other part of the trap. 

C. Transport induced ion loss 

The third step in our transport efficiency characterisation is to check that all ions leav¬ 
ing zone 1 are trapped in zone 2 by the end of the transport protocol. As only a single 
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fluorescence collecting optics is used in the experiment, the precise characterisation of the 
transfer efficiency requires that the ions are transferred back to their original position for a 
comparison between the cloud sizes. To circumvent this drawback and be able to estimate 
the one-way transfer efficiency, we identihed a transfer protocol that is efficient enough to 
serve as a standard operation. This is the case for the transfer of 100 /rs duration. Like 
mentioned previously, the transfer efficiency depends very little on the ion number as long 
as this number is larger than two thousand ions. This efficiency was estimated from several 
consecutive transport protocols to be of the order of 90% for 100 /rs. By using always the 
same protocol for the zone 2 to zone 1 transfer, we can observe the dependence of the zone 1 
to zone 2 transfer efficiency as a function of its duration, like shown on hgure 5. The results 



FIG. 5. Blue squares : Proportion of ions transferred from zone 2 to zone 1 by a transfer protocol 
of duration 100 ps and then transferred back to zone 2 by a transfer protocol of duration tg. Red 
dots (deduced from Fig. 3(b)): Proportion of ions not leaving zone 1 after transfer protocol of 
duration tg. The DC voltages are 1200 V (the lines are a guide to the eye) 


show oscillations with the transfer duration, out of phase with the proportion of ions not 
leaving trapping zone 1, which is also reproduced from Fig.3(b) on the same hgure. The 
largest two-way transfer efficiency is as high as 90% and is observed for a transport made of 
two consecutive 100 ps transport protocols. Increasing the duration of a transport protocol 
does not result in a higher transfer efficiency. Numerical simulations detailed in [22] show 
how the cloud spreading makes long transport inappropriate for large clouds. The extra 
information brought by the comparison of the two curves of Figure 5 is that the ion number 
budget evidences transfer-induced ion loss. In the following we quantify the transfer-induced 
cloud heating to look for possible correlations with the ion loss. 
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IV. EXCITATION OF MOTION 


The transport-induced motional excitation has a signature on the time evolution of the 
fluorescence. All laser frequencies are kept constant during the experiment and the temporal 
evolution of the fluorescence directly after the transport depends on the Doppler effect, which 
depends on the ion velocity along the trap axis. Heating may occur during transport, in 
that case, a re-cooling phase can be observed. The time Tf required for the fluorescence 
rate to reach its stationary value after a transfer operation is plotted in hgure 6. The chosen 
transport durations are large enough to make sure that a non-negligible proportion of the 
ion cloud arrives in zone 2. The re-cooling time Tf can vary from short (2 s) to long (10 s) 
times which shows that the Doppler shift induced by the transport depends on the duration 
of the transport. Also on hgure 6 (b) is plotted the maximum huorescence rate, which can be 
considered as a crude indication for the ion number. The graph conhrms that a fast recovery 
of the signal is not due to a lower number of ions. As seen on the comparison of the two 
curves of hgure 6, the amplitude variations of the signal are anti-correlated with the recovery 
time of the huorescence rate. In a hand-waving argument, we can interpret the time Tf as 
an indicator for the motional excitation, and deduce from this hgure that a larger number of 
ions is efficiently transferred to trapping zone 2 when this excitation is low. When a smaller 
RF amplitude is used (which results in = 0.12 instead of 0.15), the signal recovery time 
can be as small as 200 ms, showing that the ion velocities are less modihed by the transport 
process. The largest observed signal recovery time is 5 s for this lower potential value. The 
global increase of the signal recovery time with the RF-amplitude can be interpreted by 
non-linear terms in the equation of motion coupling the motion along the radial and axial 
directions, and giving rise to RF heating of the motion. In our experimental context, the 
axial potential is deformed and the non-harmonic contributions are non-negligible when the 
potential minimum crosses the site of the central electrode. The anharmonic contributions 
induce a coupling between the center of mass motion and the motion in the center of mass 
frame which is responsible for an increase of the kinetic energy of the center of mass. MD 
simulations of the transport of an ion cloud [22] showed that this contribution increases with 
the transport duration, as the cloud spreads further out. The experimental results do not 
show such a behaviour for the time-scale explored, leaving the cause of duration-dependent 
cloud heating unexplained. 
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V. TRANSFER EFFICIENCY VERSUS ION NUMBER 


We analyse the ion number effect for the transport duration which gives the highest 
two-way transport efficiency. In the described case, the protocol uses transport functions 
of duration tg = 100/is. Figure 7 shows that for clouds of less than 5 000 ions, the round- 
trip transfer efficiency increases with shrinking cloud size and can reach unity for ensembles 
of less than 2000 ions. We assume that this size effect is due to the spatial spreading of 
the cloud. This figure also shows that for this chosen transport function, the round-trip 
efficiency for shuttling is typically higher than 80 %. These high ratios can be realised with 
ion clouds of up to 10® ions. This very fast and efficient shuttling is in particular interesting 
for experiments in frequency metrology, as for example [12]. 


VI. CONCLUSION 

We have experimentally studied the transport of ion clouds in a macroscopic rf trap for 
cloud sizes as large as several tens of thousands. This transport is efficiently controlled by 
the potential applied to the central DC-electrode splitting the trap in two zones. We have 



W 



FIG. 6. a); Time Tf it takes for a transported cloud in zone 2 to recover its maximal fluorescence 
rate after a transfer from zone 1 of duration tg. b): Value of this maximal signal. The DC initial 
voltages are 1800 V. 
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FIG. 7. Ratio of the number of ions after shuttling back and forth from trapping zone 2 with two 
identical transport protocols, versus the initial number of ions. 


used a time profile for the transport function which is designed for single ion shuttling, and 
which can result in unity transfer of an ion cloud. We have also investigated the cloud’s 
response to the potential deformation indnced by the spatial translation. Our experimental 
results and their comparison with MD simulations show that the mismatch between the ideal 
and the real DC-potential prohles explains the varying temporal response of the transport 
efficiency on the duration of the transport fnnction. We can observe oscillations between 
high and low probability for an ion-cloud transfer from one trapping zone to another. The 
transfer induced heating also shows snch oscillations with a comparable time scale. 

The observed oscillations in the transfer efficiency can be modihed by choosing different 
trapping parameters, and the minimnm heating can be lowered by using smaller RF am- 
plitndes. It is therefore possible to hnd conditions for which the transfer efficiency is high 
and the motional excitation is low for the same transfer dnration. For clonds containing less 
than 2000 ions, 100 (±1.5)%-transfers can be achieved. This is another step approaching 
our objective which is to transfer large ion clonds with 100% efficiency withont heating. Onr 
best resnlts for clonds larger than 5000 ions, are transfers of 92%. 

We can tailor the temporal response to the transfer protocol by tnning the applied DC- 
voltages. It is possible to choose a ”no-retnrn” parameter set, where ions are transferred 
with a very high probability from a hrst trap to a second trap, bnt at the same time they 
do have an extremely low probability to leave the second trap. This asymmetric protocol 
allows to implement a true accumulation process, the experiment is described in [17]. 
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